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We revisit the isovector pseudoscalar sum rule determination of 7n„ + m^, 

using families of finite energy sum rules known to be very accurately satisfied 

in the isovector vector channel. The sum rule constraints are sufficiently 

strong to allow a determination of both niu + m^ and the excited resonance 
'<> ■ 

\^ • decay constants. The corresponding Borel transformed sum rules are also very 

C^ , 

well satisfied, providing a non-trivial consistency check on the treatment of 

direct instanton contributions. We obtain [m„ + tti,^] (2 GeV) = 7.8 it 1.1 MeV 



(in the MS scheme), only marginally compatible with the most recent sum 
rule determinations, but in good agreement with recent unquenched lattice 
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I. INTRODUCTION 

Because of the Ward identity d^A^^ = {rriu + rrid) ui'j^d, a study of the correlator 

n„.(g2) = t J d'xe^^- (0|T {d,A^^,ix)d.A:J (0)) |0) (1) 

allows one, in principle, to determine rriu + rud ^- A number of such studies have been 
performed ||l]-|], the most recent (Refs. (BPR) and [Q (P98)) employing, respectively, 3- 
and 4-loop expressions for the dominant D = OPE contribution. We concentrate on the 
results of P98 (which updates BPR) in what follows. P98 quotes, for m„ + mj^ 

[m„ + rud] (2 GeV) = 9.8 ± 1.9 MeV . (2) 

Recent unquenched lattice simulations, in contrast, yield [|,|1 

[ruu + md]{2 GeV) = 6.88^;^^ MeV (CP - PACS) 

[m„ + mj(2 GeV) = 7.0 ± 0.4 MeV (QCDSF - UKQCD) , (3) 

where the errors do not reflect the uncertainty involved in using perturbative versions of the 
renormalization constants. Because the consistency of the lattice and sum rule determina- 
tions is not particularly good, we revisit the sum rule treatment of Hud- 

In this paper, we study n^^ using Borel transformed and finite energy sum rules (BSR's 
and FESR's). The BSR's have the form 

M'B[UU{M') = / dse-^/^^'pUs) (4) 

Jo 

^ dse~^/'''pUs)+ dse-^/^^'p^r{s), (5) 

JO JSQ 

where M is the Borel mass, Sq the "continuum threshold", and i3 [n"J (M^) the Borel 
transform of KdiQ^) = d^U^diQ^)/ W^f). The FESR's have the form 

-z— f dsw{s)Yiud{s) = / dsw{s)pud{s) , (6) 

2m J\s\=so Jo 



^All quark masses, here and in what follows, are in the MS scheme. 



where Sq is arbitrary and w{s) is any function analytic in the region of the contour. 

The TT contribution to pud, [Pud{s)]^ = 2/^m^5 (s — m^), with f^, = 92.4 MeV, is very 
accurately known. The decay constants of the 7r(1300) and 7r(1800), needed to describe the 
remaining contributions to pud below s ~ 4 GeV^, are not known, and need to be determined 
as part of the sum rule analysis. 

The LHS of Eq. (^) can be evaluated using the OPE, provided that M is sufficiently large 
compared to the QCD scale. The condition that sq be similarly large, though necessary, is 
not sufficient for the OPE to be employed reliably on the LHS of Eq. (P) since, except for 
extremely large sq, the OPE is expected to break down near the timelike real axis p|. For 
the isovector vector (IVV) channel, this breakdown can be seen explicitly using the very 
precise spectral data available from hadronic r decay 1^: FESR's involving w{s) = s^ with 
k = 0,1,2,3 (which fail to suppress contributions from the region of the circle |s| = Sq near 
the timelike real axis) are rather poorly satisfied at scales 2 GeV < sq < rn^ [|1^. The 
breakdown of the OPE, however, turns out to be very closely localized to the vicinity of 
the timelike axis: FESR's based on weights having even a single zero at s = Sq are very 
accurately satisfied over this whole range |T^. Thus at the scales 2 GeV^ < Sq < 4 GeV^ of 
interest to us, the supplementary constraint w{so) = must be imposed in order to obtain 
reliable FESR's. We call such FESR's "pinch-weighted" , or pFESR's. 

The OPE representation of UudiQ^) is known up to dimension D = 6, with the dominant 
D = perturbative contribution known to 4-loop order |[TT| , p!2[| . Working with n"((5^), which 
allows logarithms to be summed via the scale choice /i^ = Q^, one has [p!T] , p!2 

3 {rhu + ffidf 



KdiQ': 

KdiQ'] 



D=0 



D = 4: 



{rfiu + rfid 
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1 + ^a + 14.1793a^ + 77.3683a^ 
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where a = a{Q'^) = as{Q'^)/7f, nik = mk{Q'^), with asiQ"^) and m{Q'^) the running couphng 
and running mass at scale /x^ = Q^ in the MS scheme, ^^4 and Vtf are the RG invariant 



modifications of (aG^) and {rrisss) defined in Ref. |TT[], and pvsA in Eq. (|^) describes the 
deviation of the four-quark condensates from their vacuum saturation values. We have 
dropped D = 2 contributions, which are suppressed by two additional powers of mu^d, and 
additional D = 4 contributions proportional to [niu + mrf]^m^^. The Borel transforms of 



the above expressions are well known, and may be found in Refs. [pT] , p!2 



In scalar and pseudoscalar channels, direct instanton contributions are potentially im- 



portant, but are not incorporated in the OPE representation of n^^ |13,14|. We estimate 



their size using the instanton liquid model [|I5|. ILM contributions to the theoretical side of 
the Hud BSR are given by 

3p2(m„ + mrf)2M6 



87r2 



K,{pJM'/2) + K^{pJM'/2) 



(10) 



where pi ~ (1/0.6 GeV) is the average instanton size and Ki are the MacDonald func- 
tions. ILM contributions play only a small (few percent) role in the BSR analysis at scales 
M^ > 2 GeV^, but are important for FESR analyses. For polynomial weights, ILM FESR 



contributions follow from 16 



— f dss''[Uud{s)] 
2m J\s\=so 



ILM 



-3[m„ + mj^ f^o 
An 



dss''+^jApi^)YApi^rs] ■ (11) 



We employ the following values for OPE/ILM input: pi = 1/(0.6 GeV) |Tl,|T5l; 
as{ml) = 0.334 ±.022 §; (a.G^) = (0.07 ± 0.01) GeV^ l^; {rriu + nid) {uu) = -f^ml; 
0.7 < {ss)/{uu) = r, < 1 igig; {gqaFq) = (0.8 ± 0.2 GeV^) {qq) [0; and pvsA = 
— s> 10. The D = and 4 OPE contributions are evaluated via the contour-improvement 
prescription |T^, using the analytic solutions for as(Q^) and m{Q'^) obtained from the 4- 
loop-truncated versions of the P |^ and 7 [^ functions. 



II. POTENTIAL PROBLEMS WITH THE EXISTING SUM RULE TREATMENT 

AND THE UPDATED PFESR ANALYSIS 

The BPR and P98 analyses employ FESR's with w{s) = 1, s. The global normalization of 
the resonance contributions is fixed by assuming resonance dominance of Pud at Svr threshold 
and normalizing the tails of the resonance contributions to the known ChPT threshold 
expression. Direct instanton contributions are neglected. The relative strengths of the two 
resonance contributions are constrained by optimizing a "duality" match between OPE and 
spectral ansatz versions of the ratio of the w{s) = 1- and s-weighted FESR'sQ rriu + rrid is 
then extracted from the optimized duality matching region of the w{s) = 1-weighted FESR. 
The result of Eq. (^ corresponds to sq ~ 2 GeV . The determination is, in fact, not stable 
with respect to Sq, falling roughly linearly from 9.8 MeV at Sq — 2 GeV^ to 7.5 MeV at 
So ~ 4 GeV^ (see Fig. 2 of P98). 

Two potential problems with this analysis are (1) the use of non-pinched- weighted 
FESR's at scales for which they are poorly satisfied in the IVV channel, and (2) neglect 
of direct instanton contributions. In addition, the P98 result, [mu + rridKl GeV)/[mu + 
mj(2 GeV) = 1.31 corresponds (using 4-loop running) to as{m1) = 0.307, significantly 
lower than the recent ALEPH determination. Since this will produce an overestimate of 
f^u + ^TT-d, an update of OPE input is also in order. 

Concerning the first problem, one could, of course, be lucky: the scale at which the 
OPE can be safely used right down to the timelike axis might turn out to be lower in the 
isovector pseudoscalar than in the IVV channel. If so, however, pFESR's employing the 
same spectral ansatz and same value of m„ + rrid should also be well satisfied at the scales 
used in P98. We test this possibility using pFESR's based on WN{y,A) = (1 — y)(l + Ay) 
and WD{y,A) = (1 — |/)^(1 + Ay) (where y = s/sq and A is a free parameter), which 



^The duality matching is equivalent to imposing the local duality (OPE) version of the spectral 
function, up to an overall multiplicative constant, over the whole of the matching window. 
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are known to be well satisfied in the IW channel. The resulting OPE/spectral integral 
match (corresponding to the above values for the OPE input and, as in P98, neglect of ILM 
contributions) is shown for the wn case in Fig. [1|, and is obviously quite poor. (The quality 
of the wd match is even worse.) 

We have re-analyzed n^^, using the P98 spectral ansatz as input, but fixing niu + Tn^ 
via a combined wn-, wd pFESR analysis. If we do not include ILM contributions, the 
optimized OPE/spectral integral match remains poor. Including ILM contributions produces 
a reasonable optimized match. The corresponding value of m„ + m^ is 

[ruu + mrf](2 GeV) = 6.8 MeV . (12) 

The quality of this match is shown in Fig. |^ for the wjsi family and in Fig. |^ for the wn family. 
The result of Eq. (|12]) is compatible with the P98 results corresponding to Sq — 4 GeV^ but 
significantly smaller than that corresponding to Sq = 2 GeV^. Since, in spite of optimization, 
the match for wn is best where that for wd is worst, and vice versa, it appears that some 
modification of the P98 spectral ansatz is also required. 

In this work, we aim to determine simultaneously the excited resonance decay constants, 
/^(i30o) = /i and /7r(i800) = /25 which characterize the modifications of the spectral ansatz, 
and iTiu + md. To this end, we perform a combined w^v and W£, pFESR analysis^. Our 
spectral ansatz is 

p^,{s) = 2flml5 (.s - ml) + 2flml^,,,,)B,{s) + 2flm\,,,,^B,{s) , (13) 

where Bi^2{s) are standard Breit-Wigner forms for the 7r(1300) and 7r(1800). We employ 
PDG2000 1^ values for the masses and widths. This ansatz can be used sensibly only 
up to So — ("^7r(i800) + r,r(i800)) — 4 GcV^. To maintain good convergence of the OPE 
we also require Sq > 3 GeVl For 3 GeV^ < Sq < 4 GeV^ the integrated D = OPE 



^The same treatment of the IVV channel results in a determination of /„ accurate to within the 



experimental error [22| 



series converges well for all A> 0. Larger A produces larger relative contributions from the 
resonance region, and hence aids in the extraction of /i and /2. The results of this analysis 
are 

K + ^d] (2 GeV) = 7.8 ± 0.8r ± O.^heory ± 0.4„etw MeV (14) 

f\ = 2.20 ± 0.39r ± OAStteory ± OASmethod MeV (15) 

< /a < 0.37 MeV . (16) 

The errors labelled 'T" result from varying the input resonance parameters within the 
PDG2000 errors, and are due essentially entirely to the (large) uncertainty on the 7r(1300) 
width. Those labelled "theory" reflect uncertainties in the OPE input and our estimate of 
the error associated with truncating the D = series at 0{a^). Those labelled "method" are 
obtained by studying the impact of employing different analysis windows in sq and A, and 
performing separate wn and wd analyses. Further details of the analysis, and a breakdown 
of the separate error contributions will be given elsewhere |2^. The OPE+ILM/spectral 
integral match corresponding to these results, shown in Fig. ^ for the w^ family and Fig. |^ 
for the Wd family, is obviously excellent. 

As noted above, the ILM contributions play a non-negligible role in the pFESR analysis. 
In fact, if one removes ILM contributions, an equally good OPE/spectral integral match is 
obtained, but now corresponding to [mu+mJ(2GeV) = 9.9±1.2r±1.0t/ieor,/±0.5metftod MeV, 
/i = 2.41 ± 0.50r ± 0.21theory ± 0.27„etw MeV and h = 1.36 ± 0.16r ± O.OQjwy ± 
OAlmethod McV. The pFESR analysis alone thus provides no evidence either for or against 
including ILM contributions. Fortunately, the requirement of consistency between BSR and 
pFESR analyses places non-trivial constraints on the ILM representation. This works as 
follows. The pFESR analysis provides a determination of rriu + rrid and /i^2 which is sen- 
sitive to whether or not ILM contributions are included. The output /i 2, together with 
/tt, determine the low-s part of pud, and hence can be used as input to a BSR analysis. 
The high-s part is, as usual, approximated by the continuum ansatz, with the continuum 
threshold, sq, determined by optimizing stability of the BSR output (in this case tjIu + rrid) 

8 



with respect to M. The BSR and pFESR output values for m„ + m^ should be compatible 
if the ILM representation is reasonable. Errors associated with uncertainties in input OPE 
and resonance parameter values are common to the pFESR and BSR analyses and strongly 
correlated. Additional errors are present for the BSR analysis as a result of the crudeness of 
the continuum approximation and the uncertainties in the criterion for fixing Sq. We assign 
a 20% error to continuum spectral contributions, and allow Sq to vary by ±0.5 GeV^ about 
the optimal stability value. We work in a window of Borel masses 2 GeV^ < M^ < 3 GeV^ 
for which continuum contributions are < 50% of the dominant D = OPE contribution 
and convergence of the Borel transformed D = OPE series is still good. For the case 
that the ILM contributions are included in the pFESR analysis we obtain, quoting only the 
additional errors present for the BSR analysis, 

[rriu + nid] (2 GeV) = 7.5 ± 0.9 MeV , (17) 



to be compared to the central value given in Eq. (|14D . The agreement is excellent. The 
stability of the BSR analysis, shown in Fig. ^ is also extremely good. In contrast, if 
ILM contributions are omitted from the pFESR analysis, the BSR result becomes [m„ + 
md]{2 GeV) = 8.8 ±0.6 MeV, incompatible with the pFESR determination. If one performs 
a pFESR optimization of /i_2 separately for each m„ ± md, one finds that, with no ILM 
contributions, the pFESR and BSR values remain inconsistent, within the additional BSR 
errors, unless the pFESR input, [niu + md]{2 GeV), is < 8.1 MeV. The corresponding 
optimized value for /i for this marginal case turns out to be consistent within errors with 
that quoted in Eq. (|15]), though the OPE/spectral integral match is significantly worse than 
that obtained for the optimized fit, including ILM contributions. The low value for mu + niij, 
thus appears to be an unavoidable feature of the combined analysis of 11^^. 

III. SUMMARY AND DISCUSSION 

Combining the 0.3 MeV difference of pFESR and BSR central values in quadrature with 
all other sources of error, we obtain, for our final result, 

9 



[m„ + md]{2 MeV) = 7.8 ± 1.1 MeV . (18) 

This is compatible, within errors, with the unquenched lattice determinations of Refs. [§|,§1, 
and with the result obtained by combining the ChPT determination R = 2ins/[mu + m^] = 
24.4 ±1.5 |^5[ with recent determinations of mg using hadronic r data [^. A recent 



summary [^] gives 83 MeV < ?7is(2 GeV) < 130 MeV, which corresponds to 6.8 < [tti^ + 
?7irf](2 GeV) < 10.7 MeV. Our analysis, in fact, favors values of nig in the lower part of this 
range, in good agreement with recent unquenched lattice results for mg. 
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FIGURES 




So (GeV=) 

FIG. 1. The wn OPE/spectral integral match corresponding to central values of all OPE 
input, the quoted P98 value [mu + md]{l GeV) = 12.8 MeV and the P98 spectral ansatz. The solid 
(dashed) lines represent the spectral (OPE) integrals. The lower, middle and upper lines for each 
case correspond to A = 0,2 and 4, respectively. 
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FIG. 2. The wn OPE+ILM/spectral integral match corresponding to central values of all OPE 
input, the ILM estimate of direct instanton contributions, and use of the P98 spectral ansatz. 
The OPE+ILM curves employ the optimized value, [mu + W'd](2 GeV) = 6.8 MeV, obtained in a 
combined w^, wd pFESR fit. The conventions for identifying spectral and OPE+ILM integrals, 
and the A = 0,2 and 4 cases, are as for Fig.|| above. 



o 
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s„ (GeV^) 

FIG. 3. The OPE+ILM/spectral integral match, as in Fig.g, except for the wd rather than wn 
weight family. 
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s„ (GeV^) 

FIG. 4. The optimized OPE+ILM/spectral integral match for the wn pFESR family, with 
niu + rriii, /i and /2 given by the central values of Eqs. (\14^, (P^) and (p!q). The labelling of the 
hadronic integrals, OPE integrals and the A = 0,2 and 4 cases, is as for Fig.|l] above. 




3.0 3.2 3.4 3.6 3.8 4.0 



So (GeV^) 

FIG. 5. The optimized OPE+ILM/spectral integral match, as in Fig. ^, except for the wd 
rather than wn weight family. 
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M" (GeV") 

FIG. 6. [mu + W'd](2 GeV), as a function of M^ for the BSR analysis described in the text. 
The soUd line corresponds to sq = 3.7 GeV , which produces optimal stability for rriu + m^ with 
respect to M^ in the window 2 GeV < M^ < 3 GeV . The lower (short) dashed line corresponds 
to So = 4.2 GeV^ and the upper (long) dashed line to sq = 3.2 GeV^. Note the compressed vertical 
scale. 
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